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Abstract 

A theorem from control theory relating the Lie algebra generated 
by vector fields on a manifold to the controllability of the dynami- 
cal system is shown to apply to Holonomic Quantum Computation. 
Conditions for deriving the holonomy algebra are presented by tak- 
ing covariant derivatives of the curvature associated to a non-Abelian 
gauge connection. When applied to the Optical Holonomic Computer, 
these conditions determine that the holonomy group of the two-qubit 
interaction model contains SU(2) x SU(2). In particular, a universal 
two-qubit logic gate is attainable for this model. 

PACS numbers: 03.67.Lx 03.65.Vf 

1 Introduction 

Controlling quantum dynamics to effect a desired unitary evolution is a fun- 
damental issue in quantum computation. Full control over the system dy- 
namics and hence the ability to realize any logic gate is called universal 
quantum computation. In recent years, there has been considerable interest 
in conditions for universality and it has been proven that for an n-level quan- 
tum system universality is a dense condition, being satisfied by almost all 
computational models 0-0. 

Despite the richness of the mathematical model, the effects of quantum 
noise make it a tremendous challenge to manifest this property in nature. 
Zanardi and Rasetti @] have proposed a novel methodology for the control 
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of quantum information which may provide a resolution to these compet- 
ing phenomena. Holonomic Quantum Computation (HQC), as introduced in 
0, is a theoretically appealing model that can provide universal computa- 
tion, and, due to the geometrical nature of the framework, possesses intrinsic 
robustness against decoherence and control imperfection. 

In this paper, we are primarily interested in the existence of the logic gates 
available to the experimentalist within the HQC framework. Specifically, we 
report a result from geometric control theory that simplifies the calculation 
of the holonomy group associated to a non-Abelian gauge connection. The 
application of this theorem to HQC establishes conditions for universality in 
general and in particular proves the universality of the Optical Holonomic 
Computer ||-[@]. We refer the reader to the literature for a more detailed 
exposition of the HQC set-up J7|, ||, techniques for the calculation of the 
holonomies [R EL and its intrinsic fault-tolerance |TT], |TT |. 



The methodology is briefly described as follows. The quantum code is 
realized by the n-dimensional eigenspace, C, of an ra-fold degenerate Hamil- 
tonian H with eigenvalue E . Let M be a ci-dimensional real parameter 
space. We suppose that the experimentalist can implement unitary trans- 
formations U[r]) depending continuously on the parameter 77 G M. A set of 
isospectral Hamiltonians is formed by the adjoint orbit of Hq, 

O{H ) = Mr,) Horftn) , v e M}. (1) 

Let 7 : [0, T] — ► M be a closed curve in parameter space. If we traverse 
this loop sufficiently slowly, the adiabatic theorem ensures that no energy 
level crossing will occur and 0{Hq) forms a family of Hamiltonians that 
drive the dynamics. Let \ip)i G C be the initial state of the system, then after 
completing an adiabatic loop in parameter space, the initial and final states 
are related by 

|^) / = e^ T r A ( 7 )|^) ? : G C (2) 

where e lE ° T is the dynamical phase which will be omitted in the following by 
setting E = 0. The matrix, 1^4(7) G U(n), is the holonomy associated to 
the loop 7 and for n > 1, these non-Abelian holonomies are the logic gates 
of the quantum computer. The holonomy or geometric phase [T2j depends 
only on the geometry of the loop 7 and can be expressed as 

r A ( 7 ) =Pexp I A (3) 
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where P denotes path ordering and the skew-symmetric matrix valued one- 
form A is known as the Wilczek-Zee connection [|13| with matrix elements 



A^^(g\U\ V )—U( V )\ g ). (4) 

In this way, Ta may be considered to be a map from the loop space of M to 
the matrix Lie group U(n). The set Hol(A):= {1^(7) | 7 G M} forms a group 
under composition of loops in parameter space and is, in general, a subgroup 
of U(n). Hol(A) is said to be the holonomy group, and the corresponding 
Lie algebra is known as the holonomy algebra. When Hol(v4) = U(n), the 
connection A is irreducible [Q. Clearly, irreducibility of the connection is 
sufficient for universal computation since any n-level unitary transformation 
may be applied to the code \ip) G C. 



2 Universality 



Since the work of Montgomery |14|, [15|, [T6|, [T?| , mathematicians and engineers 



have cast certain problems in the modelling and control of dynamical sys- 
tems subject to nonholonomic constraints in the language of gauge theory. 
These constraints come in two main varieties. A cat in free fall experiences 
a dynamical constraint set by the requirement that its angular momentum 
remain constant throughout its descent. An upside down cat with initial 
angular momentum zero endeavors to achieve a rotation and land upright 
by altering its shape. When a mechanical system interacts with its environ- 
ment, kinematic constraints are often encountered. For example, kinematic 
constraints are in force for a mobile robot with two independently controlled 
rear wheels subject to a no-slip constraint against the rolling surface ||18|| . 

Mathematically, the scenario is described by a connection on a principal 
G-bundle |19[. We recall these constructions and provide a geometric setting 



for the remarks made in the introduction. A principal G-bundle is formed by 
manifolds Q (total space), M (base space) a free Lie group action $ : GxQ — > 
Q, and the canonical projection Q — > Q/G = M. If U is a neighborhood of 
M, then Q is locally diffeomorphic to the product U x G. A smooth map 
a : U — >Q, such that 71 o a = idyi is called a local section over U. 

The fiber 7r _1 (j9) over a point p 6 Mis identical to the group orbit and 
is denoted G p . Let q denote the Lie algebra of G. For any element, (eg 
the group action $ e xp(t£) Q defines a curve though q G Q . The infinitesimal 
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generator £ g of the group action is defined as the tangent vector 



:= T+ $ exp(t£) Q- 



(5) 



t=o 



The vertical subspace V q Q is defined to be the subspace of T q Q that is tangent 
to the fiber G p , by the previous definition we have the identification V q Q = g. 

A connection A on Q is an Ad-equivariant Lie algebra-valued one-form 
A : TQ — > g such that A(£ g ) = £ for all £ G g. The horizontal space -f/^Q is 
the linear space H q Q := {X q G T 9 Q|^4 9 (X 9 ) = 0}. The local connection form 
A is defined with respect to a local section A = a* (A) . These definitions 
provide the splitting 



of the tangent vectors into horizontal and vertical components. Note that 



The projection map at a point defines an isomorphism from the horizontal 
space to the tangent space to the base space by T q 7i : H q — > T n ^M. Thus 
a curve q(t) G Q defines a curve in the base space by specifying a tangent 
vector at each point n(q(t)) = p(t). The properties of the connection and the 
uniqueness theory of ODE's provide the reverse procedure of reconstructing 
a curve in the total space given a curve in the base space called the horizontal 
lift ||19|| . Denote the horizontal lift of X G TM by X h and the horizontal 
part of Z G TQ as hZ. The horizontal lift of any closed curve in the base 
space maps the fiber to itself and corresponds to a group element g by the 
automorphism qf = & g qi- Assuming direct control over the base velocities, we 
seek a closed curve in the base space that achieves a desired group translation 
in the fiber. 

In the case of the cat (robot), the conservation law (no-slip constraint) 
defines a connection on a principal bundle with group action 5*0(3) x Q — ► 
Q (SE(2) x Q — ► Q). The key point in the modelling and control of these 
systems is that the horizontal distribution, defined by the connection, encodes 
the constraint information. Thus a curve q(t) G Q satisfies the constraints 
if its tangent vector X q (t) lies in H q ( t )Q for all t. Such a curve is called 
horizontal. Given an initial configuration % the feasibility of reaching a 
final configuration q/ is then equivalent to the existence of a horizontal curve 



T q Q = V q Q © H q Q 



(6) 



V q Q = KerT q 7T and H q Q = Ker^4 g . 



(7) 
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joining qi and qj. For control systems of this type, it is natural to define the 
reachable set from as the set of points q G Q that lie on a horizontal curve 
originating at g;. 

Theorem (Chow): Suppose Q is connected. Let X^ , 1 < i < d be a local 
frame of the horizontal space at q. Then any two points of Q can be joined by 
a horizontal curve if the iterated Lie brackets [X^ , [X^ , . . . , [X* , X£\ . . .] 
evaluated at q G Q span the tangent space T q Q for all q. 



For control systems without drift, this classical theorem [[2(J gives a suffi- 
cient condition to determine whether the reachable set is the entire manifold 
Q. If the Lie brackets defined in the theorem fail to span all of T q Q, the 
reachable set may be characterized as follows. Denote the subspace of T q Q 
defined in the theorem as 

A q = span{[X£ , [X^, . . . , [X*, X*] . . . ] ; 1 < i k < d, 1 < k < oo} (8) 

The sub-bundle A = M A q forms, by construction, an involutive distribution 
on the manifold Q. If the rank of A q is constant as q is varied, the Frobenious 
theorem [2T], |22| then asserts the existence of an integral submanifold Q C 



Q with A as its tangent space. This submanifold is invariant under the 
constrained dynamics and forms the reachable set. 

These arguments and the general principle embodied in the theorem are 
well known in the quantum computation literature. In the usual dynamical 
approach to quantum computing, the experimentalist has a repertoire of 
Hamiltonians {if;}[ =1 that act on the quantum state. If the Lie algebra 
generated by the Hi under commutation is equal to su(n) (or u(n)), then 
the system is deemed capable of performing universal computation. This 
is equivalent to the notion of complete controllability for quantum systems 
p3| , ^4] , p5| . It is not surprising, then, that Chow's theorem is decisive the 
holonomic framework as well. 

To establish universality of HQC, we must show that the holonomy group 
Hol(y4) is rich enough to generate a universal set of logic gates. The holonomy 
group is determined by the 0-valued curvature two-form defined by 

F(X ll X 2 )=dA(hX 1 ,hX 2 ) (9) 

where d denotes exterior differentiation. To evaluate the curvature we employ 
the structure equation 

^(X 1 ,X 2 ) = dA(X 1 ,X 2 ) + [A(X 1 ),A(X 2 )] (10) 
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for Xi,X 2 G T q Q. Since the connection A evaluates to zero on horizontal 
vectors, the structure equation implies that ^(X^X^) = — -4.([X^, X\ ]). 
Recalling that a connection can be defined as a projection of a vector X G 
T q Q onto the vertical subspace V q Q = g, we see that T{X h ', Y h ) is the vertical 
component of the vector [X^X^]. The Ambrose-Singer theorem expresses 
the holonomy group associated with the connection in terms of the curvature. 
As stated in p|: 

Theorem (Ambrose-Singer) : Let Q be a principal G-bundle over a manifold 
M. The Lie algebra f) of the holonomy group Hol go (^4) of a point q G 
Q agrees with the subalgebra of q spanned by the elements of the form 
J z q {X h ) Y ) where X h , Y h G H q Q and q is a point on the same horizontal 
lift as go- 

This theorem has been quoted by other authors to provide sufficient con- 
ditions for universality of HQC. The statement in italics, however, demands 
that we evaluate the curvature on the horizontal space at every point q that 
is reachable from q via a horizontal curve. This set of points, however, is 
the reachable set as defined above. 

More tractable conditions are obtained from Chow's theorem. By the 
above reasoning, elements of the form 

F(X?,X*) = -A([X?,X}]), r(x*,xi) = -A([xt,xft), ... (11) 

contribute a set of group directions obtained from brackets of horizontal 
vectors. According to Chow's theorem we must compute all the iterated 
Lie brackets of horizontal vectors. The vertical component of the vector 
[X^, [Xg, -Xg]] is given by D X hJ 7 (X2 , X3 ) and higher order Lie brackets are 
expressed as higher order covariant derivatives of the curvature. 

Corollary: Suppose Q is connected. The holonomy algebra at a point q G Q 
is spanned of by the curvature forms T{X^ x , X^ ) and the covariant derivatives 
D X h D X h . . . D X h TOQ , X£ ) evaluated at q . 



This result appears in Montgomery [14]. A proof can be found in [17 
(See also |ZE|| .) It is interesting to note that Montgomery's original moti- 
vation, in addition to the cat's problem, was the optimal control of spin 
systems. 

To apply this result to HQC we identify the relevant manifolds and the 
horizontal direction. Following Fujii || 0, let 7i be a separable Hilbert space 
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and define the manifolds 



St n (H) := {V = («!, . . . ,v n ) ETC x ... xH \ V ] V = Id nx „} 
Gr n (H) := {X G 5(H) | X 2 = X , X f = X , trX = n} 

where B{7i) denotes the set of bounded linear operators on H. These 
manifolds are known as the Stiefel and Grassmann manifolds respectively. 
They form a principal bundle with the (right) U(n) action on St n (7^) and 
the projection n : St n (7^) — > Gr n (H) given by ir(V) = VV^ . Denote 
this £/(n)-bundle by P n . Let M be the parameter space and let the map 
II : M — > Gr n (H) be given. The principal bundle of interest is then formed 
by the pullback of P by IT, Q = U*P with total space 

Q = {(A, V) e Mx St n {H) | n(A) = n(V)} 

over the base manifold M. To be precise, the left action of the matrix acting 
on a vector e C n as defined by (@) takes place in the C n vector bundle 
associated to Q |§, 0] . 

An important special case of this construction, known as the CP n model, 
has been shown to be generically irreducible In this case, 7i = C n+1 = 
and H has an n-dimensional degenerate subspace. The parameter 
space M = CP" is isomorphic to the orbit of H , 

° {Ho) ~ U(n) x U(l) ~ U{n) ~ ° P ' 

Thus II is a surjective map of M — CP 71 onto O(H ) = Gr l n+1 . Due 
to the large parameter space, this model can be shown to be irreducible 
by considering the span of the curvature form only (and not its covariant 
derivatives). Note that this model requires control over In = diniRCP™ 
parameters to control an n-level system. 

In any case, the Wilczek-Zee (|]) connection with its built in Hermitian 
structure defines the horizontal subspace by identifying horizontal vectors as 
those which are orthogonal to the fiber [13|, |J. When applied to HQC, this 
result represents a significant reduction in the control resources necessary 
for universality and thus broadens the class of quantum evolutions that are 
capable of computation. Indeed, if one considers the span of the curvature 
form only, then one incorrectly concludes that a necessary condition for 
universality of an n-level system is given by d(d — l)/2 > n 2 where d = 
dim R M m. 
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3 Optical Holonomic Computer 



It is widely believed that coherent superposition alone cannot account for 
the exponential speed-up sought in the realization of a quantum computer. 



Quantum entanglement must also be present f27fl . As observed in j|, || the 
CP™ model does not possess a multi-partite structure necessary for encoding 
entangled states. Attention is therefore directed to physical systems that 
have a multi-partite structure built in from the start and over which the 
experimentalist can exert control. While the results presented here apply to 
any HQC set-up, we are interested in a promising model coming from quan- 
tum optics where displacing and squeezing devices realize control operations 
acting on laser beams in a non-linear Kerr medium 

Let a ' , a be creation and annihilation operators of the harmonic oscillator 
and let n = a^a be the number operator. Let Tt be the Fock space generated 
by and a with basis {\u) : v = , 1 . . . }. Each qubit is encoded in the 
degenerate subspace of the interaction Hamiltonian 

H 1 = Xhn(n - 1) (12) 

where X is a constant 0. This computing scheme scales to a system of 
m qubits by employing m lasers to form the product basis \v\V2 ■ ■ ■ v m ) = 
l^i) <8> 1^2) <S> • • • I ^m) where z/j e {0, 1} . In accordance with the quantum 
circuit model, a control strategy is devised to implement all single qubit 
rotations and a non-trivial two-qubit transformation. A fundamental result 
then asserts that universality of the entire quantum register of m qubits can 
be achieved by this set of two-level local transformations [^] . 

Single mode squeezing and displacing operators are employed to control 
the single qubit, 

S(fj,) = exp (/la^ 2 — /la 2 ) D(X) = exp (Aa^ — Aa) (13) 

where fi , A G C . These operators define a two parameter orbit of Hi under 
the unitary transformation U(X, fi) = D(X)S(fi), 

OiH 1 ) = U(X, ii)H l U\\, p) . (14) 

The holonomy group associated to loops in the (A, //) parameter space is U (2) 
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To prove universality of the computational model it suffices to generate 
non-trivial U(4) holonomies. For the two-qubit system, the Hamiltonian is 
given by, 



H = Xhmim - 1 ) + xhn 2 (n 2 - V 



(15) 



where is the number operator for the z-th beam. Two-mode squeezing and 
displacing operators realize control operations, 



M(Q = exp (C a \ a 2 ~~ C a i fl 2) ^(0 = ex P (£ a i a 2 — £aia 2 



(16) 



where ( = ^e* 6 * 2 , £ 



r 3 e 



w 3 



6 C . In the adiabatic limit, the adjoint orbit 



under the action W(f , C) = N(£)M(C), 

0{H 12 )=U{^QH 12 U\^Q 



(17) 



drives the dynamics. The degenerate subspace of if 12 is given by the com- 
putational basis {|00), |01), |10), |11)}. Set |vac) = (|00), |01), |10), |11)) G 

We characterize the reachable set from qo = (m, |vac)) G Q by applying 
the conditions obtained from Chow's theorem. The local connection coef- 
ficients A v are written in terms of the base variables only, (r 2 , 2 , r 3 , # 3 ) G 
MCC 2 §,§, 



2 1 , 
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The non-zero local curvature forms F, 
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span su{2) x -u(l) . 

The block structure of these matrices suggests that new group directions 
may be obtained by taking covariant derivatives of F T2 g 2 along the base co- 
ordinate vectors ^ and -£q , 
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These matrices and the independent contributions of the previous set span 
the Lie algebra su{2) x su{2) x u(l) C u(A). The connection is not ir- 
reducible, however, the product structure of the subgroup shows that non- 
trivial U(4) transformations are attainable. This result reconciles conflicting 
results from the literature. One concludes after consideration of the cur- 
vature form only the holonomy group to be SU(2) x U(l) M. However, a 
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variant of the square root of SWAP gate, 
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can been explicitly constructed 
ment of 517(2) x SU(2). 



This universal transformation is an ele- 



4 Concluding Remarks 

Holonomic Quantum Computation represents a novel approach to quantum 
computing by employing non-Abelian geometric phases to perform informa- 
tion processing. The geometric phase is a beautiful phenomena with a long 
history in physics, mathematics, and engineering. We have shown that a 
result from control theory provides key insight into the foundations of HQC 
and opens up the possibility that more physical systems will be amenable to 
this approach. For a particular manifestation of HQC - the Optical Holo- 
nomic Computer - the results presented here provide new understanding of 
the controlled interactions attainable in the experimental set-up. Moreover, 
we have proven the existence of a universal set of logic gates. 

Acknowledgements : I would like to thank Fabio Celani, Claudio DePersis, 
and Paolo Rinaldi for helpful discussions and Professor Richard Montgomery 
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